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ABSTRACT 

Recent observations of high redshift quasars at z ^ 6 have finaUy revealed com- 
plete Gunn-Peterson absorption. However, this at best constrains the volume-weighted 
and mass-weighted neutral fractions to be x^j > 10^'^ and > 10^^ respectively; 
stronger constraints are not possible due to the high optical depth for hydrogen Ly- 
man transitions. Here I suggest certain metal lines as tracers of the hydrogen neutral 
fraction. These lines should cause unsaturated absorption when the IGM is almost 
fully neutral, if it is polluted to metallicities Z 10~'^'^ — 10~^'^Zq. Such a minimal 
level of metal pollution is inevitable in the middle to late stages of reionization unless 
quasars rather than stars are the dominant source of ionizing photons. The 01 line at 
1302 A is particularly promising: the 01 and H ionization potentials are almost iden- 
tical, and 01 should be in very tight charge exchange equilibrium with H. The SiH 
1260 A transition might also be observable. At high redshift, overdense regions are the 
first to be polluted to high metallicity but the last to remain permanently ionized, due 
to the short recombination times. Such regions should produce a fluctuating 01 and 
Sill forest which, if observed, would indicate large quantities of neutral hydrogen. The 
01 forest may already be detectable in the SDSS z = 6.28 quasar. If seen in future 
high-redshift quasars, the 01 and Sill forests will probe the topology of reionization 
and metal pollution in the early universe. If in additional the HI optical depth can be 
measured from the damping wing of a high-redshift gamma-ray burst, they will yield 
a very robust measure of the metallicity of the high-redshift universe. 



1 INTRODUCTION 

, Recent spectroscopic observations (Becker et al 2001; Pen- 
tericci et al 2001; Djorgovski et al 2001) of high-redshift 
quasars discovered by the Sloan Digital Sky Survey (Fan 
et al 2000, 2001a) have revealed long gaps in the spec- 
tra consistent with zero transmitted flux. This long-awaited 
detection of the Gunn-Peterson effect may herald the ob- 
servational discovery of the reionization epoch. However, 
the high oscillator strength of the hydrogen Lya transition 
means that complete Gunn-Peterson absorption is expected 
even for a highly ionized intergalactic medium (IGM). The 
strongest constraint comes from the Ly/3 absorption trough, 
due the weaker (by ~ 5) oscillator strength of Ly/3. From 
the z = 6.28 quasar observed by Becker et al (2001), Fan 
et al (2001b) conclude that at z ~ 6, the lower limits on 
the mass-weighted and volume-weighted neutral hydrogen 
fraction are Shi > 10~^ and Shi > 10~^ respectively, larger 
by almost two orders of magnitude from z ~ 4. Studies in- 
terpreting the observations conclude that the observed ab- 
sorption troughs are consistent with the tail end of reioniza- 
tion, or post-overlap phase after individual HII regions have 
merged (Barkana 2001; Fan et al 2001b). However, due to 
the rapid or phase-change like nature of reionization in stan- 
dard scenarios (Gnedin 2000; Razoumov et al 2001), the 



"dark ages" or pre-overlap phase, when a substantial frac- 
tion of the hydrogen in the universe was neutral, is likely not 
far off. The spectra of the z — 6.28 quasar suggests a very 
rapid evolution in the effective optical depth and thus the 
ionizing radiation field and effective neutral fraction (Fan et 
al 2001b). This implies that a slightly higher redshift quasar 
may indeed lie within the pre-overlap era. 

Unfortunately, even if such a quasar is discovered, we 
may not learn anything new about the pre-reionization 
epoch. The hydrogen Lyman-series absorption trough sat- 
urates fully for a neutral hydrogen fraction at mean den- 
sity xm ~ 10~*; because the transmitted flux declines ex- 
ponentially with an increasing neutral fraction, we do not 
have the power to distinguish between an almost fully neu- 
tral IGM and one with only a tiny neutral fraction. As the 
IGM becomes almost fully neutral we might observe the 
red damping wing of the Gunn-Peterson trough (Miralda- 
Escude 1998). Unfortunately, the highly luminous quasars 
presently observed probably ionize their surroundings on 
several Mpc scales; the consequent reduction in optical 
depth precludes observation of the red damping wing (Cen 
& Haiman 2000; Madau & Rees 2000). The only hope of 
detecting a damping wing would be to discover objects that 
ionize only a small region of the surrounding IGM: either 
a high-redshift gamma-ray burst (which has a very short 
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duty cycle) or less luminous quasars or galaxies (which can 
be detected through gravitational lensing, Ellis et al 2001). 
Alternatively, one might hope to detect the Lya halo sur- 
rounding a high-redshift source as Lya photons scatter and 
redshift in the surrounding neutral IGM (Loeb & Rybicki 
1999). This also suffers from the difficulty that sources tend 
to ionize their surroundings; furthermore, the low surface 
brightness of the halo implies that detection is likely only 
possible with NGST. 

What can be done with present-day technology? 
Clearly, we need absorption-line probes which are still un- 
saturated when the IGM is predominantly neutral. This is 
possible if the absorbers are much less abundant than hy- 
drogen or have very small oscillator strengths. They should 
have ionization potentials similar to that of hydrogen in 
order to trace the HI fraction as faithfully as possible. In 
addition, their absorption lines must lie redward of the hy- 
drogen Lya wavelength A = 1216 A, in order to avoid con- 
fusion with the lower-redshift Lya forest. In this paper, I 
suggest metal absorption lines as a probe of the neutral 
IGM. Metals are a natural probe: for a fully neutral IGM, 
THhya ~ 10^ at z ~ 6, aud while the oscillator strengths of 
metal UV/optical transitions (/ ~ 10~^ — 1) are roughly 
comparable to that of hydrogen Lya, the abundance by 
number of metals should be lower by ~ 10~® — 10~^, im- 
plying Tinotais ~ lO"'^ — 1. The most uncertain aspect of 
this calculation is the degree to which an IGM polluted by 
metals can still remain neutral. I argue that because over- 
dense regions are the first to be polluted with metals but 
the last to be permanently ionized (due to the short recom- 
bination time), a scenario of a neutral but metal-polluted 
IGM is plausible. Nonetheless, because of this uncertainty, 
a null detection of absorption will only yield a constraint on 
the joint metallicity/ionization state of the IGM. A positive 
detection, however, may be our best hope of unveiling an 
almost fully neutral IGM with observations of high-redshift 
quasars in the near future. In all numerical estimates, I as- 
sume a ACDM cosmology with {^Im ,^A,^b, h, agf^-i ,n) = 
(0.35, 0.65, 0.04, 0.65, 0.87, 0.96). 



2 CAN METALS BE SEEN AT HIGH 
REDSHIFT? 

It is useful to begin by ruling out some promising possibili- 
ties. The best absorption line probes would involve promor- 
dial elements, which are not afflicted with uncertainties asso- 
ciated with the (unknown) high-redshift metal abundance. 
The most obvious candidate, hydrogen 21 cm absorption, 
has too weak an oscillator strength, the optical depth across 
a Hubble volume is r = 4.2 x 10"^ (a;Hi) (^^p) {^^Y^^ 
(where Ts is the spin temperature), and the observational 
difflculties in detecting a signal are formidable^ (Shaver et al 
1999). Similarly, H2, which lacks a dipole moment, has too 
low an oscillator strength. HD does have a dipole moment 



However, 21 cm emission from the neutral IGM might be de- 
tectable with the Square Kilometer Array (Tozzi et al 2001) 



and higher oscillator strengths (by a factor of ~ 1000) but in- 
sufficient to offset its low primordial abundance xhd ~ 10~^. 
Likewise, the optical depth of lithium is appreciable only at 
high redshift z ~ 500 (Loeb 2001). Deuterium has roughly 
the right abundance (~ 10~^) and oscillator strength (same 
as H), but its Lya wavelength lies too close to the H Lya 
transition (offset by only ^ 82km s~^) to be useful. Hel has a 
meta-stable state 2^S state , which becomes populated dur- 
ing the recombination cascade. Resonance-line absorption 
from this state occurs at long wavelengths ~ 4471, 5876 A, 
which are redward of hydrogen Lya as required. However, 
due to the relatively short lifetime, ~ lO^s, of this state, it 
is appreciably populated only in highly dense and signifi- 
cantly ionized gas; n(2^S)/n(He+)= 5.8 x 10"''T4"^-^V(1 + 
3110T4-°'^n-^) « 1.9xlO~''r4~° V (Clegg 1987) , where rie 
is the electron number density in cm"''. The only possibility 
left is metal lines. 

The optical depth of a line across a uniform IGM which 
has been homogeneously polluted by metals is: 

.^0.16.,( ^° ,) f^) f^) C-±±Y'' (1) 

V2.7 X 10-v Vo.o5y V1302A/ V 7 J ^' 

where Xi is the fraction of the metal atoms a at the ap- 
propriate ionization state i, Xa = Z/Zq x (ria/nn)© is the 
abundance by number of metal a relative to hydrogen, / 
and A are the oscillator strength and rest wavelength of the 
appropriate transition. I adopt (nc,©, ^0,0, nsi,©, npo,©) = 
(3.58, 8.49, 0.33, 0.295) x lO"'*nH,0 for the solar abundance 
of carbon, oxygen, silicon and iron respectively (Anders & 
Grevesse 1989). For some metals, particularly Si, this may 
be an underestimate: the supernovae of supermassive stars 
which are thought to form out of very low/zero metallicity 
gas overproduce a elements such as Si, S and Ca by factors 
of a few compared to solar ratios (Heger & Woosley (2001); 
see Fig. 1 of Oh et al (2001)). In Table |l[ I list various lines 
which I have identified as promising tracers of the IGM ion- 
ization state^: those with ionization potentials close to that 
of hydrogen and strong resonance lines redward of hydrogen 
Lya. If the IGM is polluted to a metallicity of Z ~ W^ '^Zq 
at z ~ 6, the optical depths are fairly high: there would be a 
ffux decrement of ~ 5 — 20% blueward of these lines, which 
is certainly detectable in high signal-to-noise spectra. Note 
that ions with ionization potential li < 13.6eV may be rare 
as the universe is optically thin to radiation at these wave- 
lengths and these atoms are very easily ionized to the next 
stage. Submillimeter fine structure lines of metals have os- 
cillator strengths which are too low for absorption to be 
detectable. Although (as we shall see) metal line absorption 
likely produces a fiuctuating forest rather than a mean fiux 
decrement, the relative values of r provide a good estimate 
of the importance of various transitions. 

01 is a particularly promising tracer. Its ionization po- 
tential /, = 13.618eV is only AE = 0.19eV higher than 
that of hydrogen : therefore a detection of 01 almost cer- 
tainly signals the presence of neutral hydrogen. In fact, oxy- 
gen should be locked in tight charge exchange equilibrium 

^ Tables for the atomic constants are available at 
http: / /www.pa.uky.cdu~vcrncr/lines.html 
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Table 1. Metal absorption lines which may potentially be ob- 
servable in a nearly neutral IGM. is the ionization potential 
of the ion, A and / are the absorption wavelength and oscillator 
strength, and t is the optical depth of the line across a uniform 
IGM at 2 = 6.3, assuming a metaUicity of Z = IO'^-^Zq. Only 
the species with !{ > 13.6eV (OI, Fell, Sill) are likely to be abun- 
dant, since the universe is optically thin to radiation below the 
Lyman limit. 



Ion 


h (eV) 


A(A) 


/ 


t{z = 6) 


Fel 


7.87 


2484 


0.557 


0.11 


Sil 


8.1 


2515 


0.236 


0.05 


CI 


11.2 


1657 


0.148 


0.23 


OI 


13.6 


1302 


0.05 


0.14 


Fell 


16.1 


2383 


0.3 


0.05 


Sill 


16.34 


1260 


1.18 


0.13 



with hydrogen, through the processes O + — > + H°, 
0+ H° -> O -f H+ (Osterbrock 1989). The equilibration 
timescale is ~ ~ 1-7 x IO^xhiA (^-7^) years 

(where A is the gas overdensity) , much shorter than the 
Hubble time. Therefore the OI fraction should be very 
accurately given by — | ^^^exp (AiJ/fcaT), where 

exp {AE/ksT) 1 for'/csT > AS = 0.19eV. 

The Sill 1260 A line is another promising absorption 
feature. Note that Sill has a 1304 A transition at almost the 
same wavelength as the 01 transition, but due to the weak 
oscillator strength / = 0.0871 it has an optical depth only 
7% that of 01, unless the Si to O ratio is strongly enhanced. 
Similarly, the A = 1526 A transition of Sill (/ = 0.132) has 
r « O.llroi. 

There are two large uncertainties in the above esti- 
mates. The first is the mean metaUicity of the IGM at 
high redshift. The very short timescale on which massive 
stars evolve ~ 10^ — 10*^ years implies that the IGM could 
have been polluted very early. The metal abundance of the 
quasar environment inferred from spectra of the z — 6.28 
SDSS quasar is indistinguishable from that of lower redshift 
quasars; there is little or no evolution in the observed su- 
persolar metallicities from z ~ 6 to 2 2, implying that the 
first stars around quasars must have formed at z > 8 (Pen- 
tericci et al 2001). The mean metaUicity of the IGM at 2 « 3 
at the lowest observable column densities is Z ^ 10~^'^Zq 
(Songaila 1997; Ellison et al 2000). From the spectra of 32 
quasars in the redshift range 2.31 — 5.86, Songaila (2001) 
finds no evolution in the mean universal metaUicity of the 
IGM in the redshift range z = 1.5 — 5.5. The minimum 
value she finds at 2 = 5 is Z > 10~^'^Zq; this is a strictly 
minimal estimate because it assumes CIV and SilV are the 
dominant ionization stages (which may no longer be true at 
higher redshifts as the metagalactic radiation field drops and 
gas densities increase); furthermore, lines are substantially 
undercounted at the highest redshifts because of high noise 
levels at the longest wavelengths. Very plausible theoretical 
scenarios can be constructed in which an early generation of 
stars pollute the intergalactic medium to mean metallicities 
Z ^ 10~^-^ Zq with volume filling factors > 20% by 2 ~ 9, 
without significant hydrodynamic perturbation of the IGM 
(Madau, Ferrara & Rees 2001). 

In fact, we emphasize that for reionization to take place. 



a significant amount of star formation and thus metal ejec- 
tion must take place. Madau & ShuU (1996) find that the 
energy in the ionizing continuum released by stars is ~ 0.2% 
of the rest-mass energy of metals produced, or ~ 1.8 MeV 
per metal baryon; this is fairly independent of the IMF since 
the same massive stars which produce ionizing photons pro- 
duce metals. A similar relation holds for supermassive stars 
M > IOOMq thought to form out of metal-free gas, which 
explode as pair-instability supernovae: they eject ~half of 
their mass as metals (Heger & Woosley 2001), and re- 
lease 1.8 X lO''* HI ionizing photons s ^ Mg over a life- 
time ~ 3 X lO'^yrs (Bromm, Kudritzki & Loeb 2001), which 
yields ~ 3.3 MeV per metal baryon. Using the Madau & 
ShuU (1996) relation and assuming ~ 20 eV per HI ionizing 
photon, we can write: 

--»-<Tf£3j)(5t)+°?" 

where is the number of ionizing photons per baryon in 
the universe, Z is the mean metaUicity of the universe, /esc is 
the escape fraction of ionizing photons from their host halos, 
and nqso is the contribution from quasars (which produce 
ionizing photons but no metals). The escape fraction /esc 
is highly uncertain, with estimates ranging from ~ 5% in 
the local universe (Leitherer et al 1995; Dove, ShuU & Fer- 
rara 2000) to as high as ~ 50% in highly luminous Lyman 
break galaxies (Steidel, Pettini & Adelberger 2001). How- 
ever, even for fesc ~ 100%, the mean metaUicity should be 
reasonably high toward the tail end of the reionization pro- 
cess, Z ~ 10~^'^Zq. Another uncertainty is the filling fac- 
tor of metal-polluted regions fz ■ We show in §^ that most 
reasonable values of fz should give rise to an absorption 
signal. Significant retention of metals by their host halos is 
unlikely due to the shallow potential wells predominant at 
these early epochs; in addition metal-enriched material is 
much more easily ejected from halos than the ambient gas 
(MacLow & Ferrara 1999). The only case where the universe 
can be reionized without significant co-production of metals 
is if quasars are the dominant ionizing source. Again this 
is highly uncertain, but note that the comoving emissivity 
of quasars at 2 = 5, as inferred from the quasar luminos- 
ity function, is insufficient to keep the universe reionized at 
2 = 5 by an order of magnitude (Madau, Haardt & Rees 
1999; Fan et al 2001a); the density of quasars at high red- 
shift is also constrained by the lack of faint red unresolved 
objects in the Hubble Deep Field (Haiman, Madau & Loeb 
1999). Except for the very early stages of reionization (which 
are unlikely to be accessible with present-day instruments, 
in any case), the IGM is likely to be polluted to sufficiently 
high metaUicity to make metal-line absorption studies feasi- 
ble. 

The second, much larger uncertainty is whether regions 
which are pre-enriched with metals can still remain neu- 
tral. As previously noted, the filling factor of ionized re- 
gions should be considerable once the IGM is polluted up 
to metallicities Z ^ 10~^'^ — W~^'^Zq, unless the escape 
fraction is very small /esc < 1%, and it is possible that 
all metal-polluted regions will also be ionized. Indeed, for 
fesc > few %, the typical size of HII regions at 2; ~ 9 will 
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be greater than that of the metal-laden supernovae-driven 
superbubble (Madau, Ferrara & Rees 2001); an ionization 
front precedes the metal-pollution front. Even if the ionizing- 
photon escape fraction is extremely small, the metal pollu- 
tion front is likely to coUisionally ionize the IGM by shock 
heating it to T > 10'''^K, due to the high speed of the ex- 
panding superbubble. 

However, it is important to realize that hydro- 
gen recombination times at high redshift, tree « 3 x 



1000 



the Hubble time, 



V 10 J 



>0.7 



V 10* K 

ionization fraction, Xe ~ 1/(1-1- {t/trac)) 



yrs, are short compared to 



-1.5 



V 0.7 



SO the 



0.3A"^ (txt) " (toTk) ' ^^"^ could become 

70% neutral. The lifetime of sources is likely to be short: 
the lifetime of massive stars is Ims ~ 10® — lO^yrs, and 
the duty cycle of quasars is probably of order the Edding- 
ton timescale ~ lO^yrs (such a lifetime is consistent with 
current observations of quasars; see Blandford (1999) and 
references therein). Early reionization in the pre-overlap era 
is likely a highly stochastic process in which regions of the 
IGM are ionized, polluted with metals, and then recom- 
bine and become largely neutral until another source lights 
up. While early reionization is temporary and the ionization 
state of any given region of the IGM fluctuates, metal pollu- 
tion is a permanent process and the metallicity of the IGM 
rises monotonically. Note also that although high-density fil- 
aments may initially be coUisionally ionized by the accretion 
shock during gravitational collapse, the gas will eventually 
cool and recombine. 

In particular, overdense regions are the first to be pol- 
luted up to high metallities (due to their proximity to sites 
of star formation) but they are the last to remain perma- 
nently ionized (due to the short recombination times). From 
numerical simulations, Cen & Ostriker (1999) find that 
metallicity depends very strongly on local density: at every 
epoch, higher-density regions have much higher metallicities 
than lower-density regions. In fact, the highest-density re- 
gions quickly saturate at near-solar metallicities early on. 
These results are in much better agreement with observa- 
tions than scenarios in which metal pollution is uniform. 
Miralda-Escude, Haehnelt & Rees (2000) point out that 
in an inhomogeneous universe reionization should begin in 
voids and gradually penetrate into overdense regions; the 
regions of highest density are the last to be reionized. This 
picture is strongly substantiated in numerical simulations of 
reionization (Gnedin 2000). These arguments suggest that 
a line of sight to a high-redshift quasar will intersect regions 
at or above the mean density which are largely neutral but 
nonetheless polluted with metals. Such overdense regions are 
the most likely sites to produce the metal absorption lines 
we seek. 

The temperature of metal-polluted gas is likely to be 
~ 10* - lO^K (Madau, Ferrara & Rees 2001). It can be 
shock heated up to ~ lO^K by the expanding superbub- 
ble, but cools rapidly by Compton cooling ofi'the CMB on a 



timescale tcomp = 2.3 x 10* {^) yrs. Below T - lO^'K the 



gas recombines and Compton cooling (as well as hydrogen 
line cooling) is no longer efi'ective; the gas then cools only 




log(Z/Zj 

Figure 1. Overdensity A; up to which the universe is ionized 
in the post-overlap era, as a function of the mean metallicity 
of the universe Z, assuming a (fairly robust) relation between 
the metals and ionizing photons produced by massive stars. Sohd 
lines are for z = 6, dashed lines for z = 9. The curves are plotted 
for different values of /esc/(l — /qso) where /osc is the escape 
fraction of ionizing photons from star-forming halos, and /qso is 
the fractional contribution of quasars to the ionizing background. 
The relation between A; and Z shown here is used in Figure [sj 



by adiabatic expansion on the Hubble expansion timescale. 
In highly overdense and metal-polluted regions metal-lines 
cooling will become important, but these correspond to col- 
lapsed halos which are in any case unstable to star forma- 
tion. T ~ 10'*K is also the equilibrium temperature if a 
photoionizing background is present. 



3 THE OI FOREST 

3.1 A simple model for gas clumping and metal 
pollution 

Let us now quantify the effects of gas clumping in the IGM. 
Miralda-Escude, Haehnelt & Rees (2000) find the following 
to be a good fit to the probability distribution by volume of 
gas overdensities A seen in the LCDM numerical simulations 
of Miralda-Escude et al (1996): 



Pv(A)dA = Aexp 



(A- 



-2/3 



Co) 



2(25o/3)2 



A-^dA 



(3) 



where they tabulate values for A,l3,Co,5o at different red- 
shifts z=2, 3, 4 and 6. One can extrapolate their results to 
higher redshifts by using So = 7.61/(1 -I- z) (which fits their 
results to better than 1%), assuming /3 = 2.5 (correspond- 
ing to an isothermal slope for high-density halos), and fixing 
A, Co by requiring the total mass and volume to be normal- 
ized to unity. Their fit is valid if the gas is smoothed on the 
Jeans scale for a gas temperature T ~ lO^K; we have argued 
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that high-redshift metal-polluted gas should indeed be at ap- 
proximately this temperature. The fraction of baryons above 
a given overdensity Ai by volume and by mass is then given 
by /v(A,) = /^Pv(A)dA and /m(AO = APv(A)dA 
respectively. 

In order to compute the optical depth of the IGM 
to metal line absorption we need to specify two unknown 
functions, the metallicity Z{A, z) and ionization fraction 
Xi{A, z) of the IGM as a function of overdensity and redshift, 
but only in the combination Y{A, z) = Xi{A, z)Z{A, z). We 
can make progress by making some simplifying assumptions. 
The growth of metallicity has two free parameters: the mean 
metallicity of the universe Z (a proxy for the total amount 
of star formation), and the volume filling factor of metal- 
poUuted regions fz- These two parameters are obviously 
inter-related, but due to the large uncertainties we treat 
them as independent free parameters, with an upper bound 
on the filling fraction fz < 0.2(Z/10~^-^Zq)^''^ (where the 
normalization is based on the model of Madau, Ferrara & 
Rees (2001), and the exponent mimics the energy depen- 
dence of the adiabatic Taylor-Sedov solution). Filling factors 
as low as ~ 1% are possible if the metal-enriched ejecta have 
magnetic fields which resist mixing with the IGM; metal- 
polluted regions could be restricted to magnetized "streaks" 
which are then sheared and distorted by subsequent gravita- 
tional clustering (Madau, Ferrara & Rees 2001). The growth 
of reionization has two free parameters, the filling factor Q of 
ionized regions, and the overdensity up to which gas is ion- 
ized Ai. Again these parameters are obviously related, but 
we can make the approximation that they are decoupled, 
with the following argument. The early stages of reioniza- 
tion are characterized by reionization of the voids (Ai < 1) 
and growth of the filling factor of ionized regions Q. How- 
ever, at some point overlap occurs {Q ~ 1), and overdense 
regions start to be reionized, and Ai grows (this hinges on 
the fact that high-density regions only occupy a small frac- 
tion of the volume). As noted by Miralda-Escude, Haehnelt 
& Rees (2000), the value of ^°^<='^'^p at which Q ~ 1 depends 
on the nature of the ionizing sources: for dim but numerous 
sources /\°^'='^1''p ^ whereas for bright but rare sources 
^overlap larger, since higher density regions have to be 
ionized before percolation can occur. We therefore treat it 
as a free parameter. 

The progress of reionization and the growth in metallic- 
ity Z are coupled via equation (^. In particular, the relation 
between (Q, Ai) and Z depends on the escape fraction of ion- 
izing photons fesc and the relative contribution of QSOs to 
reionization /qso = n^^'^ /n-y. In the early stages of reion- 
ization, recombinations are unimportant and Q oc oc Z. 
In the late stages when recombinations dominate the con- 
sumption of ionizing photons, we can relate Z and the over- 
density up to which gas is ionized Ai by equating the num- 
ber of ionizing photons per baryon with the mean number of 
recombinations per baryon in a Hubble time of the ionized 
gas: 

n^{Z) =tHanCmiiA,) (4) 

where CHii(Ai) = JJ^' A^Pv(A)dA is the clumping factor 
of ionized gas. This assumes that most star formation and 



hence metal pollution occured during the last Hubble time. 
The relationship between Ai and Z is shown in Fig |l|, for 
different values of /gsc/(1 — /qso): the smaller the value of 
this parameter, the larger the amount of star formation and 
thus metal pollution Z needed to keep the universe ionized. 

With this simple picture we can make an ansatz for the 
evolution of Y{A,z). I assume that metal pollution begins 
in the most overdense regions, while reionization begins in 
the most underdense regions. At any given epoch, I assume 
Z Zcrit for A > A{^ and Z ^ otherwise. Similarly, 
I assume the neutral fraction xhi ~ 1 for A > Af and 
xhi ~ otherwise. The value of A(^ is given by the implicit 
equation /v(Ai) — fz, while Zcrit = Z/ fM{A{^). The value 
of Af is Af < 1 during the pre-overlap phase and is given 
by equation (^) in terms of Z, /esc, /qso in the post-overlap 
phase, when recombinations are important. There are two 
limits to consider. In the early stages of reionization, when 
Af^ < A{^ , the high-density regions where metals reside 
are largely neutral, Xi ~ 1. The growth in Y with time 
is dominated by the increase in metallicity. This epoch can 
therefore be characterized by the two parameters {Z, fz) - In 
the late (post-overlap) stages of reionization, when Af^ > 
A(^ , the evolution in Y is dominated by the evolution of 
Xi, when increasingly dense (and metal-polluted) regions be- 
come ionized. In this regime the model can be specified in 
terms of the parameters (Z, /z, /csc/(l — /qso)), and Af 
can be computed from equation (^ . The transition between 
these two regimes occurs when the volume filling fraction 
of metal polluted regions and neutral regions are compara- 
ble, /z ~ 1 — Q. Since we expect fz < 0.2 (metal pollution 
should not be effective in voids, which occupy most of the 
volume), the transition regime occurs roughly at the point 
of overlap, when Q 1. The transition occurs very quickly: 
1 — Q evolves extremely rapidly at the point of overlap, when 
the mean free path of ionizing photons rises on a very short 
timescale (of order the light travel time across an HII region) 
and the ionizing background increases dramatically. By con- 
trast, Z evolves on the timescale for structure formation or 
tn- In fact, Z may drop at the point of overlap since the 
sudden rise in the IGM temperature and Jeans mass could 
cause a drop in the comoving star-formation rate (Barkana 
& Loeb 2000). 

In summary, for most of the gas, the parameter Y = XiZ 
rises in the early stages of reionization as the metallicity 
grows, peaks at an epoch roughly correponding to the over- 
lap epoch, and then falls as the gas is reionized. Note that 
for higher A, the peak value of Y is higher and occurs at 
progressively later epochs, since gas at higher overdensities 
is only reionized at later times, and so there is a longer time 
interval for metal pollution to take place. We now examine 
the observational predictions of this simple model. 

3.2 Observational Predictions 

We begin by computing the mean Gunn-Peterson absorp- 
tion in a clumpy and inhomogeneously polluted universe. 
The optical depth due to regions of overdensity A is given 
by r{A,z) — ^ ■^(^'^) °^'^f^''V o(Z, Xi, z), where To is the op- 
tical depth of the line in a uniform IGM which is uniformly 
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Figure 2. Number of OI lines above a given column density Nqi 
observable in the pre-overlap phase, when metal-polluted high 
density regions are still largely neutral, as a function of the vol- 
ume filling factor of metals fz- Solid lines are for z = 6.5, dashed 
lines for 2 = 9. The column density scales directly with the as- 
sumed mean metallicity Z. The slope of the relation depends on 



whether A 



Noi 



< A 



fz 



(in which case Nn^^g increases with fz) 
( decreases with fz)- See text for details. 



The number of observable Sill lines with comparable equivalent 
widths is roughly half that of OI. 



polluted to a metallicity Z and has a mean ionization frac- 
tion Xi. The mean metal- line Gunn-Peterson absorption A 
due to gas with A > Acrit is then given by: 



(A) 



)P(A)dA. 



(5) 



poo 

y4(AcHt) = (l - e-") = / (1 - 

Note that in general ^( Acrit) is smaller than A for a uniform 
IGM. For instance, for Acrit ~ 3 at z = 6 and Z = 10"^ '^ ^0, 
we have Tctt ~ 0.02a;i, rather than r ~ O.lAxi for a uniform 
IGM. Tcff falls exponentially with increasing Ai, due to the 
exponentially small fraction of baryons at high overdensities. 

The tight charge-exchange equilibrium between OI and 
HI implies that there is a direct relation between their ef- 
fective optical depths, independent of gas clumping or the 
nature of the ionizing radiation field: 



eft 



= 1.1 X 10~ 



10-2 



efl 



(6) 



where (Z) is the HI column-density weighted metallicity 
of the universe (as opposed to the mean metallicity Z; in 
general {Z) > Z). Therefore, if we could measure both 
Tof and Tnf , we can obtain a robust and relatively model- 
independent measure of the metallicity of the high-redshift 
universe. Such a fortuitous occasion might arise if we could 
observe a high-redshift gamma-ray burst, which does not ex- 
hibit a strong proximity effect, and therefore allows measure- 
ment of r^f by measuring the shape of the damping wing. 



The transmitted flux recovers its full value at AA/A ^ 0.1 
redward of the damping wing (Miralda-Escude 1998), so 
a clean separation of the contribution of OI absorption at 
1302(1 + Zs)A (AA/A ~ 0.07) might be possible, particu- 
larly if the shape of the damping wing is well-constrained. 
Although 01 absorption will probably produce a fluctuat- 
ing forest, Tof can be obtained by smoothing the spectrum. 
The ability to measure a mean 01 decrement of a few percent 
depends on the accuracy to which sky lines and absorption 
from other sources can be ruled out (see discussion at end 
of section). In the absence of a measurement of rgf , a con- 
servative lower bound on the metallicity of the high-redshift 
universe can still be placed from equation (0) by assuming 
the IGM to be fully neutral, Xi ^ 1, and uniform (since 
(roi) > rSf ). 

Metals in the high-redshift IGM will probably not pro- 
duce a Gunn-Peterson-like absorption trough but rather a 

forest of metal lines which particularly in the case of 

01 will provide a snapshot of neutral regions along the 

line of sight. Schaye (2001) shows that many properties of 
the Lya forest can be understood by associating the charac- 
teristic lengthscale of absorbers with the local Jeans length. 
This allows us to associate a column density for an ion i 
with a given overdensity A: 



;.2 X 10^^ cm^ 



(2.7 X 10-6 ) 



1 + z 
7 

1/2 / A 



1/2 



(7) 



/Ay/2 

UO^K/ V3/ 



.lO^K, 

where Xi and Xa are the ionization fraction and metal num- 
ber abundance. This corresponds to an equivalent width 

where I have 



5.8X10- (- 



(0.05) (" 



used the fact that lines will always be on the linear portion 
of the curve of growth. This gives an observed equivalent 
width 



(8) 



which is certainly detectable with extended integration on 
Keck. At a given overdensity A the OI 1302 A and Sill 
I26OA equivalent widths are roughly equal; the increased 
oscillator strength of the Sill line compensates for its re- 
duced abundance. Since on average Wi increases monoton- 
ically with A, for any given Wi, there exists some Ai such 
that A > Ai ^ W > Wi. The mean spacing between lines 
with > HA can be estimated by the mean separation be- 
tween contours of overdensity Ai in the universe. This can 
be estimated as (Miralda-Escude, Haehnelt & Rees 2000): 



A, = A„ [1- Fv{A,)]~^^'\ 



(9) 



where Fv{Ai) is the fraction of the volume with A < Ai, 
and XoH — 60kms-^ (basically determined by the Jeans 
length) is a good fit to numerical simulations. 

There is only a finite stretch of the spectrum over which 
01 and Sill absorption can be seen before it becomes con- 
fused with the hydrogen Lya forest. Suppose we observe 
a bright quasar which ionizes its surroundings so that the 
damping wing of the Gunn-Peterson trough does not ex- 
tend redward of the Lya line. A photon can redshift for 
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I = 5X/X ~ 20, 000 km s ^ from the OI absorption edge 
at 1302(1 + Zs)A and I ^ 10,000kms"^ from the Sill ab- 
sorption edge at 1260(1 -|- Zs)A, before it encounters the HI 
Gunn-Peterson trough (this interval is actually somewhat 
smaller, due to the finite width of Lya and NV emission 
lines). The average number of lines with N > Ni is ^ Xi/l, 
while the probability that no lines with N > Ni are seen is 
exp{—\i/l). 

In Figure ^ we show the mean number of OI lines above 
a given column density A^oi detectable in the pre-overlap 
phase. In this phase, {Z, fz) are free parameters. For a given 
fz, the column densities can be rescaled to the assumed 

mean metallicity via A'^oi = .'Voi ^ j^Q-if ^ ■ We see that 

~few lines can be seen in the lO^** — lO^^cm"^ range in the 
pre-overlap era for filling factors fz ~ 10%. Somewhat more 
lines can be seen at higher redshift since A^'i oc (1 -I- z)^^^ . 
The two distinct slopes in the relation can be easily under- 
stood. For a given Z, increasing the metal-filling factor fz 
increases the number of patches along a line of sight to a 
quasar which are metal-polluted, but decreases their mean 
metallicity Zcrit. For a low column density threshold, the 
former effect dominates and Mines increases with fz; for a 
high column density threshold, increasing fz increases the 
number of lines which fall below threshold, and Nuncs de- 
creases with fz- 

Similarly, in Fig. |^ we plot the number of lines which 
can be seen in the post-overlap era at redshift 2 = 6. The re- 
lation between Z and Ai shown in Fig. |l| has been assumed. 
The different lines illustrate the effect of varying the model 
parameters around the fidicial model (./ese, /z, A^oi'') = 
(0.1, 0.1, lO^^cm"^). Again, the change in slope can be easily 
understood: as the ionized overdensity Ai increases the num- 
ber of lines seen initially increases, because of the larger im- 
plied star-formation rate and thus higher Z. At some point 
the decrease in filling factor of neutral regions overwhelms 
the increase in metallicity, and the number of lines decreases. 
The overdensity Ai at which this break occurs depends on 
/esc/(l — /qso) and fz, since these parameters control the 
relationship between Ai and .Zcrit = Z/fz- For instance, 
the universe is ionized up to a much higher Ai for a given 
Z if /esc/(l — /qso) is large. We see that it is quite plau- 
sible for us to see 01 absorption lines with observed equiv- 
alent widths Wx « 5.3 A (^qT5^^^) (i±^) A in the SDSS 
z = 6.28 quasar, which lies just at the tail end of reionization 
Ai ~few. 

Interesting (though model-dependent) constraints on 
/esc , /qso and fz might be possible from measurements 
of the OI forest. For instance, the length of the dark re- 
gion and upper limit on the transmitted flux in the Lya, (3 
troughs give constraints on Ai for a given model of struc- 
ture formation; Fan et al (2001b) find for their z — 6.28 
quasar that Ai ~ 3. From Figure hi we see that this implies 
Z ^ W-^-^Zq (^) ^iz^so^ ^ £g^jj.iy jjjgjj metallicities. 

From Figure ^ the number of OI lines above a given column 
density in the post-overlap phase depends on /z,/esc; the 
observed number might constrain their value. The universe 
is likely to be in the pre-overlap era if a very long Gunn- 
Peterson trough with no detectable fiux is seen. From Fig- 
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Figure 3. Number of OI lines with Nqi > lO^^cm"^ observable 
in the post-overlap era at z = 6 as a function of Ai , the overden- 
sity to which the IGM is assumed to be ionized. The fidicial model 
(shown in bold) is for fz = 10%, /esc = 10% (for a substantial 
QSO contribution, /esc should be replaced with /esc/ (1 — /qso))- 
Other lines show the eff'ect of varying /zi/esc and Nqi- The 
break in the slope of the relation can be easily understood; see 
text for details. The number of observable Sill lines with compa- 
rable equivalent widths is roughtly half that of OI. 

ure ^ the number of detectable lines above a given column 
density might constrain Z,fz- The numerical value of the 
constraints are of course very model-dependent and should 
not be over-interpreted. Still, some interesting statements 
might still be made with reasonable confidence. For instance, 
if no OI lines can be seen in the post-overlap era despite a 
deep Gunn-Peterson damping trough, then Z is low: this 
implies the escape fraction of ionizing photons is close to 
unity and/or quasars are the dominant ionizing source. Al- 
ternatively, the majority of the metals are highly ionized, 
either because most metals reside in voids (although this 
is unlikely given the finite speed at which metal pollution 
fronts can propagate), or the high-density regions in which 
they reside are constantly illuminated by an ionizing source. 

The main observational obstacle to detecting the OI 
forest is confusion with other sources of line absorption. In- 
trinsic absorption within the quasar host can be constrained 
by the relative absence of absorption features in the rest 
frame 1216 - 1302A range of the large sample of lower red- 
shift quasars. Indeed, because of the reasonably wide wave- 
length interval in which OI and Sill absorption can be seen 
(~ 20, OOOkms"^ and ~ 10, OOOkms"^ respectively), at the 
shorter wavelengths within this range intrinsic absorption 
can be ruled out. Broad absorption line (BAL) quasars could 
have winds which show such features, though usually the 
outflows are highly ionized. Confusion with long wavelength 
metal lines such as CIV or Mgll from lower-redshift systems 
is another problem. Such an origin can be constrained by 
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the absence of a damped-Lya system at the corresponding 
redshift in the low redshift Lya forest. The third and po- 
tentially most serious problem is the fact that the night sky 
becomes increasingly noisy at these near-IR wavelengths, 
and the atmospheric OH forest becomes important. The ac- 
curacy with which the OI forest can be detected therefore 
depends on the accuracy with which telluric features can be 
divided out via a standard-star calibration. Also, there are 
stretches between the night sky OH forest lines in which no 
absorption should be seen, so OI absorption can be identi- 
fied if it falls within these regions. Ultimately, the presence 
of neutral gas can be corroborated with simultaneous de- 
tections of 01, SiH and possibly Fell absorption features. 
A larger quasar sample which shows long stretches of com- 
plete HI Gunn-Potorson absorption and in which the density 
of OI and Sill absorption features is higher in higher redshift 
quasars should be an unambiguous signature of almost fully 
neutral patches of gas at high redshift. 



4 DISCUSSION 

The SDSS 1030-1-0524 quasar at ^; = 6.28 shows tantahzing 
absorption features blueward of the 01 1305A line (Becker 
et al 2001). There also appears to be a fairly deep ab- 
sorption feature blueward of the Sill 1260A. Could the ab- 
sorption lines described in this paper been seen already? 
Unfortunately, the quasars SDSS 1044-0125(2 = 5.80), 
0836-F0054(2 = 5.82),1306-h0356(z = 5.99) also show some 
absorption features in the same wavelength interval; in par- 
ticular, SDSS 1306-1-0356 shows a very strong absorption 
feature at ~ 7130A, with no detected flux over ~ 80A 
(this has been tentatively identified as CIV absorption at 
z = 4.86). These features cannot correspond to 01 absorp- 
tion lines: unless they correspond to regions of anomalously 
high metallicity, the associated hydrogen column densities 
would be Nhi > 10^°cm~^, and all flux at the hydrogen 
Lya wavelength should be obliterated, while some flux is 
still seen there. These lines are probably associated with 
metal lines (e.g. Mgll) from lower redshift absorbers, and 
illustrate a generic difficulty in observing the features pro- 
posed in this paper. On the other hand, 01 lines are not 
ruled out in the z — 6.28 quasar because of the complete 
damping of flux at Lya, Ly/3 wavelengths. As we have seen, 
these absorption features can still arise in the post-overlap 
epoch when regions with A > few are largely neutral; in- 
deed, up to a few absorption lines with observed equiva- 
lent widths W\ ~ 5 {^^) ( Kii-i'lm-^ ) ^ might be seen. 
Note that the spectra of Becker et al (2001) have been 
smoothed to 4A pixel" ^. The absorption features blueward 
of OI 1305A cannot be definitely identified as 01 absorp- 
tion. Two of them lie at the same wavelengths as bright sky 
emission lines and are probably due to imperfect sky sub- 
traction. A third line with observed frame equivalent with 
~ 25 A is a possible candidate; however, it lies very close to 
the rest frame 01 wavelength and could be due to intrinsic 
absorption. More cannot be said without careful study of the 
spectrum. A definitive detection of the OI forest can prob- 
ably only be done with much higher signal-to-noise spectra 
of the same quasar. 



The estimates in this paper can be addressed with in 
much greater detail with numerical simulations. In particu- 
lar, I used very simple ansatzes for the dependences of metal- 
licity and ionization fraction with overdensity, Z{A),Xi{A) 
which in fact should be highly stochastic and spatially vary- 
ing. They can be much better modelled in a self-consistent 
fashion in simulations which attempt to model the metal 
pollution (Cen & Ostriker 1999; Aguirre et al 2001) and 
radiative transfer (Gnedin 2000; Razoumov et al 2001), 
particularly since the rise in metallicity and ionization frac- 
tion are inter-related. The spatial structure of the OI forest 
can be computed by shooting lines of sight through a simu- 
lation box. If the OI forest is indeed seen, such studies will 
be urgently needed to provide a more realistic interpretation 
of the observations. 

The scenario in this paper is not significantly altered 
if a substantial X-ray background due to high redshift su- 
pernovae (Oh 2001) or quasars (Venketesan et al 2001) is 
present. The large mean free path of hard photons means 
that they can ionize the IGM fairly uniformly, but beyond 
Xe ~ 0.1 most of the energy of an energetic electron created 
goes into Coulomb heating the gas rather than coUisional 
ionization (Slmll & van Steenberg 1985); predominantly 
neutral regions should therefore still exist. 

We will gain a wealth of information about early metal 
pollution and the reionization process if the 01 and Sill 
forests are seen. They will be direct probes of the topology 
and history of gas clumping, metal pollution and reioniza- 
tion in the early universe. As we have seen, if we assume a 
relation between the metals and ionizing photons produced 
by massive stars, they could potentially also provide indirect 
constraints on the escape fraction of ionizing photons from 
star forming lialos and the QSO contribution to the ionizing 
background, and the filling factor of metal pollution. They 
may also give clues as to the nature of the ionizing sources: 
the structure of forest lines should look different if the uni- 
verse were reionized by rare but luminous source as opposed 
to abundant but faint sources, since in the former case higher 
overdensity regions have to be ionized before overlap can be 
achieved. 01 and HI will be locked in very tight charge ex- 
change equilibrium at high redshift. If we are lucky enough 
to observe the rest-frame optical afterglow of a high-redshift 
gamma-ray burst and measure both rgf (from the damping 
wind) and Tof , we will have a direct measure of the mean 
metallicity of the universe at high redshift, independent of 
gas clumping or the form of the ionizing radiation field. Oth- 
erwise, a lower limit on the metallicity from measurement of 
Tof alone is possible. A null detection of the OI, Sill forests 
will yield constraints on the parameter Y'i(A) = Xi{A)Z{A), 
but a positive detection will be tremendously exciting and 
almost certainly signal the presence of almost fully neutral 
hydrogen at high redshift. To date, the OI and Sill forests 
may be our only probes of nearly neutral gas in the pre- 
reionization epoch observable with current technology. 
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